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ABSTRACT pressure and mean velocity profile data. The location of 

the boundary layer separation was dependent of the 
The paper experimentally studies the effects of Reynolds number. It is observed that starting point of 

periodic unsteady wake flow and different Reynolds the separation bubble and the re-attachment point move 

numbers on boundary layer development, separation and further downstream by increasing Reynolds number 

re-attachment along the suction surface of a low from 110,000 to 150,000. Also, the size of the 

pressure turbine blade. The experimental investigations separation bubble is smaller when compared to that for 

were performed on a large scale, subsonic unsteady Re=l 10,000. 

turbine cascade research facility at Turbomachinery 

Performance and Flow Research Laboratory (TPFL) of NOMENCLATURE 

T exas A&M University. The experiments were carried 

out at Reynolds numbers of 110,000 and 150,000 c blade chord 

(based on suction surface length and exit velocity). One c axial chord 

steady and two different unsteady inlet flow conditions q pressure coefficient, C = - ?i 

with the corresponding passing frequencies, wake d R rod diameter P (?r ” PXni 

velocities, and turbulence intensities were investigated. shape factor, H ] .=6/6, 

The reduced frequencies chosen cover the operating h m maximum separation huhhle height 

range of LP turbines. In addition to the unsteady L ss suction surface length 

boundary layer measurements, surface pressure hi number of samples 

measurements were performed. The inception, onset, N number of wake cycles 

and the extent of the separation bubble information p, static pressure taps i=l,.., 4 8 

collected from the pressure measurements were Ps>Pt static and total pressure at the inlet 

compared with the hot wire measurements. Re LSS Reynolds number based Re = N 

The re suits presented in ensemble -averaged, and the S b blade spacing 

contour plot forms help to understand the physics of the § r ro( j spacing 

separation phenomenon under periodic unsteady wake s streamwise distance from the leading edge of 

flow and different Reynolds number. It was found that the blade 

the suction surface displayed a strong separation bubble Sq streamwise distance from the leading edge to 

for these three different reduced frequencies. For each the trailing edge of the blade 

condition, the locations defining the separation bubble s md maximum separation bubble height at a 

were determined carefully analyzing and examining the 
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streamwise distance from blade leading edge 
re-attachment point of the separation bubble at 
a streamwise distance from blade leading edge 
starting point of the separation bubble at a 
streamwise distance from blade leading edge 
time 

turbulence intensity 
belt translational velocity 
axial velocity 
exit velocity 
velocity 

fluctuation velocity 

blades stagger angle 

boundary layer thickness 

boundary layer displacement thickness 

boundary layer momentum thickness 

boundary layer energy thickness 

kinematic viscosity 

cascade solidity, a = c/S B 

one wake-passing period 

flow coefficient, (p - VJU 

Zweifel coefficient 
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INTRODUCTION 


In recent years, gas turbine engine aero dynami cists 
have focused their attention on improving the efficiency 
and performance of the low pressure turbine (LPT) 
component. Research works from industry, research 
centers, and academia have shown that a reduction of 
the blade number can be achieved without substantially 
sacrificing the efficiency of the LPT blading. This 
reduction contributes to an increase in thrust/weight 
ratio, thus reducing the fuel consumption. Contrary to 
the high pressure turbine (HPT) stage group that 
operates in a relatively high Reynolds number 
environment, dependent on operation conditions, the 
LPT experiences a variation in Reynolds number 
ranging from 50,000 to 250,000. Since the major 
portion of the boundary layer, particularly along the 
suction surface is laminar, the low Reynolds number in 
conjunction with the local adverse pressure gradient 
makes it susceptible to flow separation, thus increasing 
the complexity of the LPT boundary layer 
aerodynamics. The periodic unsteady nature of the 
incoming flow associated with wakes that originate from 
upstream blades substantially influences the boundary 
layer development including the onset of the laminar 
separation, the extent of the separation bubble, and its 


turbulent re-attachment. Of particular relevance in the 
context of LPT aerodynamics is the interaction of the 
wake flow with the suction surface separation bubble. 
While the phenomenon of the unsteady boundary layer 
development and transition in the absence of the 
separation bubbles has been the subject of intensive 
research, that has led to better understanding the 
transition phenomenon, comprehending the multiple 
effects of mutually interacting parameters on the LPT 
boundary layer separation and their physics still requires 
more research. 

The significance of the unsteady flow effect on 
efficiency and performance of compressor and turbine 
stages was recognized in the early seventies by several 
researchers. Fundamental studies by Pfeil and Herbst 
[1], Pfeil et al. [2], and Orth [3] studied and quantified 
the effect of unsteady wake flow on the boundary layer 
transition along flat plates. Schobeiri and his co-workers 
[4], [5], [6], [7] experimentally investigated the effects 
of the periodic unsteady wake flow and pressure 
gradient on the boundary layer transition and heat 
transfer along the concave surface of a constant 
curvature plate. The measurements were systematically 
performed under different pressure gradients and 
unsteady wake frequencies using a squirrel cage type 
wake generator positioned upstream of the curved plate. 
Liu and Rodi [8] carried out the boundary layer and heat 
transfer measurements on a turbine cascade, which was 
installed downstream of a squirrel cage type wake 
generator mentioned previously. 

Analyzing the velocity and the turbulence structure 
of the impinging wakes and their interaction with the 
boundary layer along the concave side of the curved 
plate, Chakka and Schobeiri [7] developed an 
intermittency based unsteady boundary layer transition 
model. The analysis revealed a universal pattern for the 
relative intermittency function for all the frequencies and 
pressure gradients investigated. However, the curved 
plate investigations were not sufficient to draw any 
conclusion with regard to an eventual universal 
character of the relative intermittency function. Further 
detailed investigations of the unsteady boundary layer 
on a high Reynolds number turbine cascade by Schobeiri 
et al. [9], [10] and its subsequent analysis [11] and [12] 
verified the universal character of the relative 
intermittency function. For this purpose, Schobeiri et al. 
[9] utilized a conceptually different type wake 
generator, which is also used for the investigation 
presented in this paper. Fottner and his coworkers [13], 
[14] and Schulte and Hodson [15] used the same wake 
generating concept for the investigations on the 
influence of the unsteady wake flow on the LPT- 
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boundary layer. Kaszeta, Simon, and Ashpis [16] 
experimentally investigated the laminar-turbulent 
transition aspect within a channel with the side walls 
resembling the suction and pressure surfaces of a LPT 
blade. Lou and Hourmouziadis [17] investigated the 
mechanism of separation, transition, and re -attachment, 
and the effect of oscillating inlet flow conditions on 
laminar boundary layer separation along a flat plate 
under a strong negative pressure gradient which was 
similar to the LPT pressure gradient. This was simulated 
by contouring the top wall. They studied the Reynolds 
number effect on the transition region. Their results 
showed that the higher Reynolds numbers cause an 
earlier transition and reduction of the transition length, 
while the separation point does not change its location. 
Using the top wall contouring, Volino and Hultgren 
[1 8] performed an experimental study and measured the 
detailed velocity along a flat plate which was subjected 
to a similar pressure gradient as the suction side of a 
low pressure turbine blade. They also stated that the 
location of the boundary layer separation does not 
strongly depend on the Reynolds number or free- stream 
turbulence level, as long as the boundary layer remains 
non-turbulent before separation occurs. Furthermore, 
they showed that the extent of the transition is strongly 
dependent on the Reynolds number and turbulence 
intensity. 

Using the surface mounted hot film measurement 
technique, Fottner and his coworkers [13] and [14], 
Schroder [19], and Flaueisen, Flennecke, and Schroder 
[20] documented strong interaction between the wakes 
and the suction surface separation bubble on the LPT 
blades, both in the wind tunnel cascade tests and in a 
turbine rig. Furthermore, they investigated the 
boundary layer transition under the influence of the 
periodic wakes along the LPT surface and found that 
the interaction of the wake with the boundary layer 
greatly affects the loss generation. Shyne et al. [21] 
performed an experimental study on a simulated low 
pressure turbine. The experiments were carried out at 
Reynolds numbers of 100,000 and 250,000 with three 
levels of free-stream turbulence. They indicated that the 
transition onset and the length are strongly dependent 
on the free-stream turbulence. As the free-stream 
turbulence increases, the onset location and the length 
of the transition are decreased. Treuren et al. [22] 
performed an experimental study along a LPT surface at 
the very low Reynolds number of 25,000 and 50,000 
with different free-stream turbulence levels. They 
showed that a massive separation at the very low 
Reynolds number of 25,000 is persistent, in spite of an 
elevated free stream turbulence intensity. However, at 


the higher Reynolds number of 50,000, there was a 
strong separation on the suction side for the low free- 
stream turbulence level. The separation bubble was 
eliminated for the higher free-stream turbulence level of 
8-9%. The investigations by Halstead et al. [23] on a 
large scale LP turbine uses surface mounted hot films to 
acquire detailed information about the quasi-shear stress 
directly on the blade surface. Investigations by 
Cardamone et al. [14] and Schroder [19] indicate that 
the benefit of the wake -boundary layer interaction can 
be used for the design procedure of modem gas turbine 
engines with a reduced LPT blade number without 
altering the stage efficiency. 

Most of the studies mentioned above on LP turbine 
cascade aerodynamics have largely concentrated on the 
measurement of the signals stemming from hot films 
mounted on the suction and pressure surfaces of the 
blades under investigation. Although this technique is 
qualitatively reflecting the interaction of the unsteady 
wake with the boundary layer, because of the lack of an 
appropriate calibration method, it is not capable of 
quantifying the surface properties such as the wall shear 
stress. The few boundary layer measurements are not 
comprehensive enough to provide any conclusive 
evidence for interpretation of the boundary layer 
transition and separation processes and their direct 
impact on profile loss, which is a critical parameter for 
blade design. Furthermore, the numerical simulation of 
the unsteady LPT blade aerodynamics using 
conventional turbulence and transition models fails if it 
is applied to low Reynolds number cases. Recent work 
presented by Cardamone et al. [14] shows that in the 
steady state case at Re= 60,000, the separation is 
captured, however, for the unsteady case, the separation 
bubble is not reproduced. 

A recent experimental study by Schobeiri and 
Oztiirk [24], [25] investigated the physics of the 
inception, onset and extent of the separation bubble 
along a low pressure turbine blade which was the first 
part of a series of investigations carried out at TPFL. A 
detailed experimental study on the behavior of the 
separation bubble on the suction surface of a highly 
loaded LPT blade under a periodic unsteady wake flow 
was presented in [24]. Surface pressure measurements 
were performed at Re= 50,000, 75,000, 100,000, 
125,000. Increasing the Reynolds number has resulted 
in no major changes to the surface pressure distribution. 
They concluded that the unsteady wake flow with its 
highly turbulent vortical core over the separation region 
caused a periodic contraction and expansion of the 
separation bubble. It was proposed that, in conjunction 
with the pressure gradient and periodic wakes, the 
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temporal gradient of the turbulence fluctuation, or more 
precisely the fluctuation acceleration d V rm /dt 
provides a higher momentum and energy transfer into 
the boundary layer, energizing the separation bubble and 
causing it to partially or entirely disappear. They found 
that f° r dv rms /dt > 0 , the separation bubble starts to 
contract whereas for dv rms /dt < 0 , it gradually 
assumes the shape before the contraction. They argued 
that not only the existence of higher turbulence 
fluctuations expressed in terms of higher turbulence 
intensity influences the flow separation, but also its 
gradient is of crucial importance in suppressing or 
preventing the onset and the extent of the separation 
bubble. They stated that the fluctuation gradient is an 
inherent feature of the incoming periodic wake flow and 
does not exist in a statistically steady flow that might 
have a high turbulence intensity. They also stated that, 
unsteady wake flow with its highly turbulent vortical 
core over the separation region, caused a periodic 
contraction and extension of the separation bubble and 
a reduction of the separation bubble height. Increasing 
the passing frequency associated with a higher 
turbulence intensity further reduced the separation 
bubble height [25]. 

The objective of the present study dealing with the 
specific issues of the LPT boundary layer aerodynamics 
is to provide a detailed unsteady boundary flow 
information to understand the underlying physics of the 
inception, onset, and extension of the separation bubble 
for different Reynolds numbers. Furthermore, the 
unsteady boundary layer data from the present and 
planned experimental investigations will serve to extend 
the intermittency unsteady boundary layer transition 
model developed by Schobeiri and his co workers [7,11, 
12] to the LPT cases, where a massive separation occurs 
on the suction surface at a low Reynolds number at the 
design and off-design incidence. Furthermore, the 
experimental results are intended to serve as benchmark 
data for a comparison with numerical computation using 
DNS or RANS-codes. 

It is well known that the boundary layer 
measurement is one of the most time consuming 
aerodynamic measurements. Any attempt to increase the 
number of parameters to be studied would inevitably 
result in substantial increase of the measurement time. 
Considering this fact, the research facility described in 
[9] and [10] with state-of-the-art instrumentation has 
been substantially modified to study systematically and 
efficiently the influence of the periodic unsteady and 
highly turbulent flow on the LPT cascade aerodynamics 
at the design and off-design incidence angles, where the 
Reynolds number, wake impingement frequency, free- 


stream turbulence, and the blade solidity can be varied 
independently. 

EXPERIMENTAL RESEARCH FACILITY 

To investigate the effect of unsteady wake flow on 
turbine and compressor cascade aerodynamics, 
particularly on unsteady boundary layer transition, a 
multi-purpose large-scale, subsonic research facility was 
designed and has been taken into operation since 1993. 
Since the facility in its original configuration is described 
in [9], [10] and [24] only a brief description of the 
modifications and the main components is given below. 
The research facility consists of a large centrifugal 
compressor, a diffuser, a settling chamber, a nozzle, an 
unsteady wake generator, and a turbine cascade test 
section as shown in Figure 1 . The compressor with a 
volumetric flow rate of 15 m 3 /s is capable of generating 
a maximum mean velocity of 100 m/s at the test section 
inlet. The settling chamber consists of five screens and 
one honeycomb flow straightener to control the 
uniformity of the flow. 

Two-dimensional periodic unsteady inlet flow is 
simulated by the translational motion of an unsteady 
wake generator (see Figure 2), with a series of 
cylindrical rods attached to two parallel operating timing 
belts driven by an electric motor. To simulate the wake 
width and spacing that stem from the trailing edge of 
rotor blades, the diameter and number of rods can be 
varied. The rod diameter, its distance from the LPT 
blade leading edge, the wake width and the 
corresponding drag coefficient is chosen according to 
the criteria outlined by Schobeiri et al. [26]. The belt- 
pulley system is driven by an electric motor and a 
frequency controller. The wake-passing frequency is 
monitored by a fiber-optic sensor. The sensor also 
serves as the triggering mechanism for data transfer and 
its initialization, which is required for ensemble- 
averaging. This type of wake generator produces clean 
2-dimensional wakes, whose turbulence structure, decay 
and development is, to a great extent, predictable [26]. 
The unsteady boundary layer transition and heat transfer 
investigations [9-12] performed on this facility serve as 
the bench mark data for validation of turbulence models, 
transition models, and general code assessments. 

T o account for a high flow deflection of the LPT 
cascade, the entire wake generator and test section unit 
including the traversing system, was modified to allow 
a precise angle adjustment of the cascade relative to the 
incoming flow. This is done by a hydraulic platform, 
which simultaneously lifts and rotates the wake 
generator and test section unit. The unit is then attached 
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Figure 1. Turbine cascade research facility with the components and the adjustable 


test section 



Figure 2. Wake Generator 
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Table Is Parameters of turbine cascade test section 


Parameters 

Values 

Parameters 

Values 

Inlet velocity 

^ = 4 m/s 

Inlet turbulence intensity 

Tu in = 1.9 % 

Rod translational speed 

U= 5.0 m/s 

Blade Re-number 

Re =110,000 

Nozzle width 

W = 200.0 mm 

Blade height 

h B = 200 mm 

Blade chord 

c = 203.44 mm 

Cascade solidity 

G= 1.248 

Blade axial chord 

c^ = 182.85mm 

Zweifel coefficient 

1|J A = 1.254 

Blade suction surface length 

L ss = 270.32 mm 

Cascade angle 

(p = 55° 

Cascade flow coefficient 

$ = 0.80 

Cascade spacing 


Inlet air angle to the cascade 

a, = 0° 

Exit air angle from the cascade 

cc 2 = 90° 

Rod diameter 

D r = 2.0 mm 

Rod distance to lead, edge 

L r = 122 mm 

Cluster 1 (no rod, steady) 

S R = 00 mm 

Q - parameter steady case 

Q = 0.0 

Cluster 2 rod spacing 


Q - parameter for cluster 1 

Q = 1.59 

Cluster 3 rod spacing 


Q - parameter for cluster 2 

Q = 3.18 


to the tunnel exit nozzle with an angular accuracy less 
than 0.05°, which is measured electronically. 

The special design of the facility and the length of 
the belts (L belt = 4,960 mm) enables a considerable 
reduction of the measurement time. For the present 
investigation, two clusters of rods with constant 
diameter of 2 mm are attached to the belts as shown in 
Figure 2. The two clusters with spacings S R = 160 mm 
and S R = 80 mm are separated by a distance which does 
not have any rods, thus simulating steady state case (S R 
= °o). Thus, it is possible to measure sequentially the 
effect of three different spacings at a single boundary 
layer point. T o clearly define the influence domain of 
each individual cluster with the other one, the clusters 
are arranged with a certain distance between each other. 
Using the triggering system mentioned above and a 
continuous data acquisition, the buffer zones between 
the data clusters can be readily identified. 

The data analysis program cuts the buffer zones and 
evaluates the data pertaining to each cluster. 
Comprehensive preliminary measurements were carried 
out to make sure that the data were exactly identical to 
those, when the entire belt length was attached with 
rods of constant spacing, which corresponded to each 


individual cluster spacing. The cascade test section 
shown in Figure 1, located downstream of the wake 
generator, includes 5 LPT blades with a height of 200.0 
mm and the chord of 203.44 mm. For boundary layer 
investigations, five identical “PakB” airfoils designed by 
Pratt & Whitney were implemented whose cascade 
geometry is given in Table 1. 

The blade geometry resembles the essential feature 
such as the laminar boundary layer separation that is 
inherent to typical LPT blades. The blade geometry was 
made available to NASA researchers and academia to 
study the specific problems of LPT flow separation, its 
passive and active control and its prevention. As shown 
in [9], this blade number is necessary and sufficient to 
secure a spatial periodicity for the cascade flow. The 
periodicity is noticed in the pressure distribution for the 
second and fourth blade for steady and unsteady flow 
conditions, and is shown in Figure 1 . These blades were 
specially manufactured for measurement of pressure and 
showed identical pressure distributions. 

A computer controlled traversing system is used to 
measure the inlet velocities and turbulence intensities, as 
well as the boundary layers on suction and pressure 
surfaces. The traversing system (see Figure 3) was 
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Figure 3. Turbine cascade research facility with 3- 
axis traversing system 


modified to allow the probe to reach all streamwise 
positions along the suction und pressure surfaces. The 
three axis traversing system is vertically mounted on the 
plexiglass side wall. Each axis is connected to a DC- 
stepper motor with an encoder and decoder. The optical 
encoder provides a continuous feedback to the stepper 
motor for accurate positioning of the probes. The 
system is capable of traversing along the suction and 
pressure surfaces in small steps up to 1 |am, and the 
third axis is capable of rotating with an angular accuracy 
less then 0.05°, which is specifically required for 
boundary layer investigations where the measurement of 
the laminar sublayer is of particular interest. 

INSTRUMENTATION, DATA ACQUISITION, 
AND DATA REDUCTION 

The data acquisition system is controlled by a 
personal computer that includes a 16 channel, 12-bit 
analog-digital (A/D) board. Time dependent velocity 
signals are obtained by using a commercial 3-channel 
(TSI), constant temperature hot-wire anemometer 
system that has a signal conditioner with a variable low 
pass filter and adjustable gain. A Prandtl probe, placed 
upstream of the diffuser, monitors the reference velocity 


at a fixed location. The pneumatic probes are connected 
to high precision differential pressure transducers for 
digital readout. Several calibrated thermocouples are 
placed downstream of the test section to constantly 
monitor the flow temperature. The wake generator 
speed and the passing frequency signals of the rods are 
transmitted by a fiber-optic trigger sensor. The passage 
signals of the rods are detected by the sensor using a 
silver-coated reflective paint on one of the belts. This 
sensor gives an accurate readout of the speed of the 
wake generator and the passing frequency of the rods. 
The signals of the pressure transducers, thermocouples, 
and trigger sensors are transmitted to the A/D board and 
are sampled by the computer. To ensure the cascade 
periodicity, the second and fourth blades are 
instrumented each with 48 static pressure taps. Two 
adjacent blades are used for boundary layer 
measurement. The taps are connected to a scanivalve, 
which sequentially transferred the pressure signals to 
one of the transducers that was connected to the A/D 
board. 

The unsteady data are taken by calibrated, custom 
designed miniature, single hot wire probes. At each 
boundary layer position, samples were taken at a rate of 
20kHz for each of 100 revolutions of the wake 
generator. The data were ensemble-averaged with 
respect to the rotational period of the wake generator. 
Before final data were taken, the number of samples per 
revolution and the total number of revolutions were 
varied to determine the optimum settings for 
convergence of the ensemble-average. 

For the steady state case, the instantaneous velocity 
components are calculated from the temperature 
compensated instantaneous voltages by using the 
calibration coefficients. The instantaneous velocity can 
be represented in the following form: 


V = V + v (1) 

where V is the mean (time-averaged) velocity and v is 
the turbulent fluctuation component. The mean velocity, 
also known as the time-average, is given by: 


V 


M 



( 2 ) 


where M is the total number of samples at one boundary 
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layer location. The root mean square value of the 
turbulent velocity fluctuation is: 



and the local turbulence intensity is defined as: 

r« toc = ^xioo = 4. (v- vf* 100 (4) 

V F N M i = 1 

For unsteady cases, the ensemble-averaged velocity, 
fluctuation velocity, and the turbulence intensity were 
calculated from the instantaneous velocity samples by: 

Hh ) * <%)> = v v (t t ) ( 5 ) 

JS j= i 

VA) - <v&)> = N ^E [ W )-<^(0>] 2 (6) 


- <Tu t (t)> = 9^ xl0 ° (7) 

where N= 1 00 is the total number of wake generator 
periods and M the number of samples taken per period. 
cVi (tj)> is the reference ensemble averaged velocity for 
the particular boundary layer traverse. 

EXPERIMENTAL RESULTS AND DISCUSSION 

T o investigate the influence of the unsteady wake 
flow on the boundary layer development along the 
suction and pressure surfaces of the LPT blade and, 
particularly, its impact on the inception and onset of the 
separation bubble, the detailed surface pressure and 
boundary layer measurements were performed at a 
Reynolds number of 110,000 and 150,000. This 
Reynolds numbers, which pertain to a typical cruise 
operation , exhibit a representative value within LPT 
operating range between 75,000 and 400,000 as 
discussed by Hourmouziadis [27]. Furthermore, it 
produces separation bubbles that can be accurately 
measured by miniature hot wire probes. For the 
Reynolds number of 110,000 and 150,000, three 
different reduced frequencies were examined. For 


generation of the unsteady wakes, cylindrical rods with 
the diameter d R = 2mm were chosen to fulfill the 
criterion that requires the generation of a drag 
coefficient C D that is approximately equal to the C D of 
the turbine blade with the chord and spacing given in 
Table 1 (for details look for the studies in [26] and 
[28]). 

To accurately account for the unsteadiness caused 
by the frequency of the individual wakes and their 
spacings, the flow velocity, and the cascade parameters, 
a reduced frequency Q is defined that includes the 
cascade solidity o, the flow coefficient (p, the blade 
spacing S B , and the rod spacing S R . Many researchers 
have used Strouhal number as the unsteady flow 
parameter, which only includes the speed of the wake 
generator and the inlet velocity. However, the currently 
defined reduced frequency Q is an extension of Strouhal 
number in the sense that it incorporates the rod spacing 
S R and the blade spacing S B in addition to the inlet 
velocity and wake generator speed. For surface 
pressure measurement rods with uniform spacings as 
specified in T able 1 were attached over the entire belt 
length. For boundary layer measurement, however, 
clusters of rods were attached, as mentioned previously. 

Surface Pressure Distributions 

Detailed pneumatic surface pressure measurements 
were taken at Re = 110,000, and 150,000. For each 
Reynolds number three different reduced frequencies, 
namely Q = 0.0, 1.59, and 3.18 are applied that 
correspond to the rod spacings S R = 80 mm, 160 mm, 
and mm. The pressure distributions in Figure 4 show 
the results of the steady case and two unsteady cases. 
The pressure signals inherently signify the time-averaged 
pressure because of the internal pneumatic damping 
effect of the connecting pipes to the transducer. The 
noticeable deviation in pressure distribution between the 
steady and unsteady cases, especially on the suction 
surface, is due to the drag forces caused by the moving 
rods. The drag forces are imposed on the main stream 
and cause momentum deficiency that lead to a reduction 
of the total and static pressure. 

The time-averaged pressure coefficients along the 
pressure and suction surfaces are plotted in Figure 4. 
The suction surface (upper portion), exhibits a strong 
negative pressure gradient. The flow accelerates at a 
relatively steep rate and reaches its maximum surface 
velocity that corresponds to the minimum C p = -4.0 at 
s/s c = 0.42. Passing through the minimum pressure, the 
fluid particles within the boundary layer encounter, a 
positive pressure gradient that causes a sharp 
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Re=1 10,000 


Re=1 50,000 




Figure 4. Static pressure distributions at two different Re-numbers and reduced frequencies Q=0, 1.59, 3.18 
(no rod, 160 mm, 80 mm), SS=Separation start, SE= Separation end 


deceleration until s/s 0 = 0.55 has been reached. This 
point signifies the beginning of the laminar boundary 
layer separation and the onset of a separation bubble. As 
seen in the subsequent boundary layer discussion, the 
separation bubble characterized by a constant C p - 
plateau extends up to s/s 0 = 0.746, thus occupying more 
than 19% of the suction surface and constituting a 
massive separation. Passing the plateau, the flow first 
experiences a second sharp deceleration indicative of a 
process of re-attachment followed by a further 
deceleration at a moderate rate. On the pressure surface, 
the flow accelerates at a very slow rate, reaches a 
minimum pressure coefficient at s/s 0 = 0.42 and 
continues to accelerate until the trailing edge has been 
reached. Unlike the suction surface, the pressure surface 
boundary layer does not encounter any adverse positive 
pressure gradient that triggers separation. However, 
close to the leading edge, a small plateau extending from 
s/s q — 0.08 to 0.16 indicates the existence of a small size 
separation bubble that might be attributed to a minor 
inlet flow incident angle. 

Considering the unsteady case with the reduced 
frequency Q = 1.59 corresponding to a rod spacing of 
S r =160 mm, Figure 4 exhibits a slight difference in the 
pressure distribution between the steady and unsteady 
cases. As mentioned above, this deviation is attributed 
to the momentum deficiency that leads to a reduction of 
the total and static pressure. For Re=l 1 0,000, the wakes 
have a reducing impact on the streamwise extent of the 
separation plateau. As seen in Figure 4 ( a), the trailing 
edge of the plateau has shifted from s/s c = 0.74 to s/s G = 
0.702. This shift reduced the streamwise extent of the 


separation plateau from 19% to 15% of the suction 
surface length which is in this particular case, 21% of 
reduction in streamwise extent of the separation. 
Increasing the reduced frequency to Q=3.18 by 
reducing the rod spacing to S R =80 mm causes a slight 
shift of the C p -distribution compared with Q=1 .59 case. 
One should bear in mind that pneumatically measured 
surface pressure distribution represents a time integral 
of the pressure events only. 

Increasing the Reynolds number to Re=l 50,000, 
has brought major changes in steady state Cp- 
distribution. The combination ofhigher Re -number with 
unsteady wakes reveals the noticeable deviation on the 
streamwise extent of the separation plateau. As seen in 
Figure 4 (b), the trailing edge of the plateau has shifted 
from s/s 0 = 0.74 to s/s 0 = 0.702 for Reynolds number of 
150,000. The combination ofhigher Reynolds number 
with high unsteady wakes introduce fluctuation kinetic 
energy into the boundary layer which tends to inhibit the 
separation tendency. C p - distribution clearly shows that 
the wake impingement with higher Reynolds number 
shortens the streamwise extent of the separation zone 
compared to the steady case. Also, the combination of 
higher Re-number with unsteady wakes reveals that the 
noticeable deviation in pressure distribution between the 
steady and unsteady cases discussed above is 
diminishing with increasing the Re-number as shown in 
Figure 4(a, b). Two counteracting factors are 
contributing to this deviation. The first factor is 
attributed to the momentum deficiency and the 
associated total pressure losses caused by moving 
wakes, as discussed above. The second factor pertains 
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to the energizing effect of the impinging wakes on the 
boundary layer. Although the impinging wakes cause 
velocity and momentum deficits, their high turbulence 
intensity vortical cores provide an intensive exchange 
and transfer of mass, momentum, and energy to the 
blade surface, thus energizing the low energetic 
boundary layer. In conjunction with the surface pressure 
distribution, the kinetic energy of the normal velocity 
fluctuation component plays a crucial role. In case of a 
low Re-number flow, the strong damping effect of the 
wall shear stress has the tendency to reduce the normal 
contribution of turbulence kinetic energy, thereby 
diminishing its surface pressure augmenting effect. 
Increasing the Reynolds number results in a decrease of 
the damping effect of the wall shear stress, allowing the 
kinetic energy of the normal velocity fluctuation 
component to increase the surface pressure, thus 
offsetting the wake deficit effects on the pressure 
distribution. This fact is clearly shown in Figure 4(a, b), 
where the pressure distributions of unsteady flow cases 
at Q=1.59 and Q=3.18 systematically approach the 
steady state cases at Re= 150,000. 

Detailed information regarding the structure of the 
separation bubble is delivered by means of a detailed 
unsteady boundary layer or surface pressure 
measurement by fast response probes, as will be 
discussed in the subsequent sections. 

Time Averaged Velocity Distributions 

Following the surface pressure investigations that mainly 
addressed the onset and extent of the separation zone 
discussed previously, comprehensive boundary layer 
measurements were performed to identify the 
streamwise and normal extent as well as the deformation 
of the separation zone under unsteady wake flow. The 
steady state case serves as the reference configuration. 

Consistent with the surface pressure distribution 
which is discussed above, the effect of the wake 
frequency on the time-averaged velocity profiles and 
fluctuation velocity distribution are presented for one 
steady and two unsteady inlet flow conditions on the 
suction surface along 3 1 streamwise locations for the 
Reynolds number of 110,000 and 41 streamwise 
locations for the Reynolds number of 150,000. After 
completing the velocity measurements, the boundary 
layer coordinates were transformed into a blade 
orthogonal coordinate system. Velocities at blade 
normal positions were obtained by interpolating their 
transformed values. The results showed almost no 
difference between the interpolated and non-interpolated 
velocity data. Experimental investigations were 


performed for three different values of Q = 0.0, 1.59, 
and 3.18. These values cover the reduced frequency 
range encountered in LPT-design and off-design 
operation conditions. 

The effect of wake frequency on time averaged 
velocity and velocity fluctuation distributions is shown 
in Figures 5 to 8 at 6 representative streamwise 
locations for Re=l 10,000 and Re= 150,000. Upstream of 
the separation bubble at s/s D = 0.52 and also at its 
immediate proximity s/s o =0.588, the velocity 
distributions inside the boundary layer experience a 
slight decrease with increasing the reduced frequency. 
Inside the separation bubble at s/s 0 =0.705, a substantial 
influence of the wake frequency is observed. The higher 
wake frequency introduces a fluctuation kinetic energy 
into the boundary layer trying to reverse the separation 
tendency. As it can be seen from the velocity 
distribution profiles wake impingement shortens the 
streamwise extent of the separation zone and decreases 
the bubble height, compared to the steady case, 
however, the onset of the separation bubble is not 
changed. This shows that the flow does not have the 
capability to suppress the separation bubble. It only 
reduces the separation bubble height. In the downstream 
of the separation bubble, where the flow is fully 
reattached, s/s o =0.951, the impact of the wake on the 
boundary layer is reduced. This effect is clearly shown 
in the velocity distribution at s/s o =0.95 1 . According to 
the previous investigations by Schobeiri et al. [10] on a 
HP-turbine cascade, an increased wake frequency causes 
turbulence fluctuations to rise inside 
and outside the boundary layer as shown in Figure 7 and 
8. However, in the LPT case with the boundary layer 
separation once the boundary layer is re-attached and 
the velocity distribution assumes a fully turbulent 
profile, no major changes are observed either in the 
velocity or in the fluctuation velocity distribution. 
Although Re=l 50,000 shows the same phenomenon for 
the velocity and the fluctuation velocity distribution, it 
is observed that starting point of the separation bubble 
and the re-attachment point move further downstream 
to s/s 0 = 0.56 and s/s c = 0.788 respectively. Also, the size 
of the separation bubble is smaller when compared to 
that for Re=l 10,000. 

Temporal Behavior of the Separation Zone Under 
Unsteady Wake Flow 

Velocity distributions on the suction surface with 
time as the parameter are plotted in Figure 9 to 12 for 
Re = 110,000 and 150,000 and Q = 1.59, 3.18. The 
nondimensional time (t/x) values are chosen so that they 
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Re=1 10,000, s/s o =0.502 


Re=1 10,000, s/s o =0.52 


Re=1 10,000, s/s o =0 .588 





Re=1 10,000, s/s o =0.705 



Re=1 10,000, s/s o =0.767 Re=1 10,000, s/s o =0.951 



Figure 5. Distribution of time-averaged velocity along the suction 
surface for steady case Q=0 (S R =°°) and unsteady cases Q=1.59 
(S r =160 mm) and Q=3.18 (S R =80 mm) at Re=l 10,000 


Re-150, 000, s/$ a =0.53 



Re=1 50,000, s/s 0 =0.56 Re= 150,000, s/s a =0.720 



Figure 6. Distribution of time-averaged velocity along the suction 
surface for steady case Q=0 (S R =°°) and unsteady cases Q=1.59 
(S r =160 mm) and Q=3.18 (S R =80 mm) at Re=l 50,000 
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Re=1 10,000, s/s o =0.502 


Re=1 10,000, s/s o =0.52 


Re= 110,000, s/s o =0.588 



Figure 7. Distribution of time-averaged fluctuation rms velocity along 
the suction surface for steady case Q=0 (S R =°°) and unsteady cases 
Q=1.59 (S r =160 mm) and Q=3.18 (S R =80 mm) at Re=l 10,000 


R©-150,000, s/s a =0.53 


Re=1 50,000, s/s 0 =0.56 


Re= 150,000, s/s 0 O-720 



V(m/s) V(m/s) V(m/s) 



Figure 8. Distribution of time-averaged fluctuation rms velocity along 
the suction surface for steady case Q=0 (S R =°°) and unsteady cases 
Q=1.59 (S r =160 mm) and Q=3.18 (S R =80 mm) at Re=150,000 
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represent the temporal states within one full period of 
wake passing. For Re = 110,000 Figures 9(a) to 9(e) 
show, the velocity distributions inside and outside the 
boundary layer at fixed s/s D - locations experience 
moderate to pronounced changes. Figure 9(a) represents 
the instantaneous velocity distribution upstream of the 
separation zone followed by Figures. 9(b,c,d,e) which 
represent the velocity distributions inside the separation 
zone. The last Figure 9(f) exhibits the instantaneous 
velocity distribution downstream of the separation zone. 
In discussing the following results, we simultaneously 
refer to the wake distribution as well as the turbulence 
fluctuation results. 

Figure 9(a) exhibits the velocity distribution on the 
suction surface at s/s 0 = 0.402. At this streamwise 
position, the laminar boundary layer is subjected to a 
strong negative pressure gradient. The boundary layer 
distributions at different (t/x) experience changes in 
magnitude that reflect the corresponding changes of the 
impinging periodic wake velocity. It is worth noting, 
that despite the injection of turbulence kinetic energy by 
the impinging wakes, no local instantaneous boundary 
layer transition occurs. This is because of the strong 
negative pressure gradient that prevents the boundary 
layer from becoming instantaneously transitional. 
Instantaneous velocity distributions inside the separation 
zone are shown in Figures 9(b,c,d,e,f). 

As a representative case, the results plotted in 
Figure 9(e) at s/s 0 = 0.674 is discussed. During the time 
interval from t/x close to 0.5 (1.5, 2.5, etc.) to about t/x 
= 0.75 (1.75, 2.75 etc.), the separation zone is exposed 
to the wake external flow which is under the influence 
of relatively lower turbulence. This flow does not have 
the capability to suppress the separation zone. Thus the 
separation region is clearly shown by the velocity 
distributions at t/x = 0.5 and t/x = 0.75. As the wake 
passes over the blade at s/s 0 = 0.674 introducing high 
turbulence kinetic energy into the boundary layer, the 
boundary layer is energized causing the separation zone 
to partially reduce or disappear. This leads to an 
instantaneous re -attachment. This time interval 
corresponds to the case where the flow is completely 
under the influence of wake and correspondingly the 
reattached velocity distribution assumes a turbulent 
profile characterized by the curves at t/x = 1.0, t/x = 
0.05, and t/x = 0.25 shown in Figure 9(e). To emphasize 
this statement, the steady state velocity distribution at 
the same streamwise position is also plotted in Figure 
9(e) using full circles. It shows clearly the separated 
nature of the boundary layer which coincides with the 
instantaneous velocity profile at t/x = 0.5. Intermediate 
times reflect the gradual change between the separation 


and re-attachment as the flow is undergoing the 
influence of the oncoming wake. Moving to the trailing 
edge of the separation zone, at s/s 0 = 0.705, Figure 9(f ) , 
a partial reduction in boundary layer thickness as the 
result of wake impingement is visible, however, the 
separation zone does not seem to disappear. 

Figure 10 (a) to (I) show the velocity distributions 
incident outside the boundary layer at fixed s/s 0 - 
locations experience moderate to pronounced changes 
at Q = 3.1 8. Similar instantaneous velocity distribution 
is observed when operating at a reduced frequency of Q 
= 3.18. Compared with Q = 1.59, a stronger 
suppression zone is noticed at Q = 3 . 1 8. When the wake 
passing frequency is increased over the blade, it is also 
increasing high turbulence kinetic energy. Therefore, 
boundary layer is energized comparatively more than the 
reduced frequency of Q = 1.59 and it causes the 
separation bubble partially reduce more. 

As it is seen in Figure 11 and 12, increasing the 
Reynolds number to 150,000 causes starting point of the 
separation bubble and the re-attachment point move 
further downstream.. It is also observed that increasing 
the Reynolds number reduces the size of the separation 
bubble. 

Behavior of Separation Bubble Under Wake Flow 

The effect of the periodic unsteady wakes on the 
onset and extent of the separation bubble is shown in 
Figure 13 to 16 for Re=l 10,000 and 150,000 and for 
two different frequencies, namely Q=1 .59 and Q=3.1 8. 
These figures display the full extent of the separation 
bubble and its behavior under a periodic wake flow 
impingement at different t/x. The wake propagation for 
Q=1.59 and Q=3.18 is analyzed, and the value of t/x 
corresponds to the point in the cycle at which the data 
acquisition system is triggered. During a rod passing 
period, the wake flow and the separation bubble 
undergo a sequence of flow states which are not 
noticeably different when the unsteady data are time- 
averaged. The temporal changes of the spatial position 
of the separation bubble, is compared with the time- 
averaged separation bubble. Starting with a reduced 
frequency of Q=1.59 at t/x=0.25, the separation bubble 
is under a full influence of the wake. The wake passing 
over the blade introduces a high turbulence kinetic 
energy into the boundary layer. The energized boundary 
layer bubble is partially reduced or disappeared 
compared to the time-averaged separation bubble size 
shown in Figure 13(a). As the wake passes, t/x=0.50, 
the height of the separation bubble reaches its maximum 
size at s/s o =0.705. The contraction starts again that 
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Figure 9. Distribution of the ensemble-averaged velocity development 
along the suction surface for different s/s c with time t/T as parameter 
for Q=1.59 (S r =160 mm) and Re=l 10,000 



Figure 10. Distribution of the ensemble-averaged velocity 


development along the suction surface for different s/s c with time t/T 
as parameter for Q=3.18 (S R =80 mm) and Re=l 10,000 
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Figure 11. Distribution of the ensemble- averaged velocity 
development along the suction surface for different s/s 0 with time t/T 
as parameter for Q=1.59 (S R =160 mm) and Re=150,000 



Figure 12. Distribution of the ensemble-averaged velocity 


development along the suction surface for different s/s c with time t/T 
as parameter for Q=3.18 (S R =80 mm) and Re=l 50,000 
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reduces the size of the separation bubble at t/x=0.75. At 
s/s o =0.705 the core region has slightly moved towards 
to the leading edge. At t/T=l, the full effect of the wake 
on the boundary layer can be seen before another wake 
appears and the bubble moves back to the original 
position. Once the wake starts to penetrate into the 
separation bubble, the turbulent spot produced in the 
wake paths causes a total suppression at some 
streamwise positions. Similar results are observed when 
operating at a reduced frequency of Q=3.18 and 
Reynolds number of 150,000. It is observed that 
increasing Reynolds number to 150,000 moves the 
separation bubble further downstream and the wake 
frequency exerts no influence on the position of the 
separation bubble. However, doubling the reduced 
frequency and increasing the Reynolds number is 
associated with the higher turbulence intensity that leads 
to stronger suppression of the separation bubble as 
shown in Figure 14 to 16. 

CONCLUSIONS 

A detailed experimental study on the behavior of 
the separation bubble on the suction surface of a highly 
loaded LPT blade under periodic unsteady wake flow 
were investigated. The measurements were carried 
utilizing a custom designed hot-wire probe. One steady 
and two different unsteady inlet wake flow conditions 
with the corresponding passing frequencies, the wake 
velocity and the turbulence intensities were investigated 
by utilizing a large-scale, subsonic research facility. 
Translational motion of the timing belt created the 
unsteady wake flow. It is found that the turbulent spots 
generated by the unsteady wake flow were effective in 
suppressing or reducing the size of the separation 
bubble. The results of the unsteady boundary layer 
measurements were presented in the ensemble-averaged, 
and contour plot forms. Surface pressure measurements 
were performed at Re= 110,000 and Re=150,000. At 
each Reynolds number, one steady, and two periodic 
unsteady inlet flow measurements were performed. 
Slight changes of the pressure distribution occurred, 
while operating at the unsteady flow conditions. 
Increasing the Reynolds number to Re=l 50,000, has 
brought major changes in steady state C p -distribution. 
The combination of higher Reynolds number with high 
unsteady wakes introduce fluctuation kinetic energy into 
the boundary layer which tends to reverse the separation 
tendency. C p -distribution clearly shows that the wake 
impingement with higher Reynolds number shortens the 
streamwise extent of the separation zone compared to 
the steady case. Detailed unsteady boundary layer 


measurement identified the onset and extension of the 
separation bubble as well as its behavior under the 
unsteady wake flow. Passing the wake flow with its 
highly turbulent vortical core over the separation region, 
caused a periodic contraction and expansion of the 
separation bubble and a reduction of the separation 
bubble height. Increasing the passing frequency 
associated with a higher turbulence intensity further 
reduced the separation bubble height. It is observed that 
starting point of the separation bubble and the re- 
attachment point move further downstream to s/s 0 = 0.56 
and s/s c = 0.788 respectively by increasing the Reynolds 
number to 150,000. Also, the size of the separation 
bubble is smaller when compared to that for 
Re=l 10,000. 

UNCERTAINTY ANALYSIS 

The Kline andMcClintock [29] uncertainty analysis 
method was used to determine the uncertainty in the 
velocity after calibration and data reduction for the 
single-wire probe. The Kline and McClintock method 
determines the uncertainty with a 95% confidence level. 
The uncertainty in the velocity for the single-wire probe 
after the data reduction is given in Table 3. As shown, 
the uncertainty in the velocity increases as the flow 
velocity decreases. This is due to the pneumatic pressure 
transducer having a large uncertainty during calibration. 


Table2: Uncertainty in velocity measurement for hot- 
wire probe. 


U (m/s) 

3 

5 

12 

coU/U (%) rcd 

5.78 

2.41 

1.40 
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(a) £1=1.59, t/x=0.25 (b) £1=1 .59, t/T=0.50 



Figure 13. Ensemble -averaged velocity contours along the suction surface for different s/s 0 with 
time t/T as parameter for Q=1.59 (S R =160 mm), Re=l 10,000 



Figure 14. Ensemble -averaged velocity contours along the suction surface for different s/s 0 with 
time t/T as parameter for Q=3.18 (S R =80 mm), Re=l 10,000 
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s/s 0 s/s 0 

Figure 15. Ensemble -averaged velocity contours along the suction surface for different s/s 0 with 
time t/T as parameter for Q=1.59 (S R =160 mm), Re=l 50,000 



Figure 16. Ensemble -averaged velocity contours along the suction surface for different s/s 0 with 
time t/T as parameter for Q=3.18 (S R =80 mm), Re=l 50,000 
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